Precision Measurement of Transition Matrix Elements via Light Shift Cancellation 



C. D. HeroldQ V. D. Vaidya, X. Li, S. L. Rolston, and J. V. Porto 
Joint Quantum Institute, University of Maryland and NIST, College Park, Maryland 20742 

M. S. Safronova 

Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716 

(Dated: August 22, 2012) 

We present a method for accurate determination of atomic transition matrix elements at the 1CP 3 
level. Measurements of the ac Stark (light) shift around "magic-zero" wavelengths, where the light 
shift vanishes, provide precise constraints on the matrix elements. We make the first measurement 
of the 5s-6p matrix elements in rubidium by measuring the light shift around the 421 nm and 
423 nm zeros with a sequence of standing wave pulses. In conjunction with existing theoretical and 
experimental data, we find 0.3236(9) eao an d 0.5230(8) eag for the 5s-6pi/ 2 an d 5s-6p 3 / 2 elements, 
respectively, an order of magnitude more accurate than the best theoretical values. This technique 
can provide needed, accurate matrix elements for many atoms, including those used in atomic clocks, 
tests of fundamental symmetries, and quantum information. 

PACS numbers: 32.70.Cs, 37.10.Jk, 67.85.Hj 



Precise knowledge of atomic transition strengths is 
important in many current areas of research, includ- 
ing the development of ultra-precise atomic clocks [H-El, 
studies of fundamental symmetries 0, [B[ and degenerate 
quantum gases Q, quantum information 0, Hj], plasma 
physics i, and astrophysics [a ITlj ]. For example, one 
of the largest contributions to the uncertainty budget of 
atomic clocks is the uncertainty in the blackbody radi- 
ation (BBR) shift 0, [ll|. The BBR shift is calculated 
from the difference in the electric-dipole polarizabilities 
between the clock states [l3| , and its accuracy is currently 
limited by uncertainty in atomic transition matrix ele- 
ments. The dynamic correction to the BBR shift can also 
be determined accurately if the relevant matrix elements 
are known [l[. Studies of fundamental symmetries, such 
as atomic parity violation, need transition matrix ele- 
ments to evaluate parity-violating amplitudes and vector 
transition polarizabilities 0, Accurate AMO theory 
can also be indispensable to the design and interpretation 
of experiments, where direct measurement of all relevant 
parameters is infeasible. More complicated atoms, such 
as Er 14], Dy 15 1, and Ho have recently become 



of interest, and development of new theoretical methods 
must be supported by the existence of high-precision ex- 
perimental benchmarks. 

Transition matrix elements can be difficult to measure 
or calculate accurately. State-of-the art theory predic- 
tions are often limited to a few percent uncertainty, and 
the presently available experimental techniques typically 
measure matrix elements accurately only for one or two 
of the lowest transitions. To date, the most accurate 
determinations of atomic transition matrix elements are 
through excited state lifetime measurements or photoas- 
sociation spectroscopy [l7j . The former is limited by un- 
certainty in branching ratios when multiple decay paths 
exist. As a result, although the 6p lifetime of rubidium 



was measured to 1%, no matrix elements were reported 
[lc| . The latter requires species with purely long-range 
molecular excited states where molecular theory is suf- 
ficiently well known to extract atomic properties. In 
principle, transition matrix elements can be determined 
from the ac polarizability a by measuring the ac Stark 
(light) shift of an atom exposed to light of known in- 
tensity. However, unlike dc Stark shift measurements, 
where the applied electric field can be determined geo- 
metrically [2j, it is difficult to accurately determine the 
optical intensity, limiting the efficacy of this approach. 
(Lifetime measurements avoid this calibration challenge 
by using a well known field, the vacuum.) For atoms 
with spin-dependent vector light shifts and two long-lived 
states, rf spectroscopy was used to accurately determine 
the ratio of two vector light shifts 1|| 2(j, which when 
combined with theory constrained matrix elements plj ]. 
New techniques are needed to improve accuracy, extend 
the range of measurable matrix elements, and provide 
benchmarks for theory. 

In this Letter, we present a widely applicable method, 
recently suggested in 22j |. for constraining matrix ele- 
ments at the 10~ 3 level through direct measurement of 
light shifts near "magic-zero" wavelengths (A ZC ro), where 
the combined light shift from all transitions cancels [23j . 
These A zor o wavelengths are distinct from the "magic" 
wavelengths used in optical clocks, where the light shift 
is identical for two states 24|. We measure small light 



shifts through diffraction off a standing wave, amplifying 
the diffracted populations by constructively interfering 
the effect of up to n p = 15 pulses, resulting in an rip en- 
hancement in the diffracted fraction. We make the first 
experimental determination of the rubidium 5s-6pi/ 2 and 
5s-6j>3/2 electric dipole matrix elements, to an accuracy 
of 0.3%, 10 times more accurate than the best theoretical 
values [22]. Implementation of our measurement tech- 



2 



1.5 



1000 
500 



-1000 



1 













5p 



418 



42Q 422 



424 




1.0 - 



0.5 - 



0.0 



-0.5 

400 500 600 700 800 

Wavelength (nm) 

FIG. 1. Calculated polarizability of Rb (atomic units of ag). 
The vertical dashed lines are at the np-5s transitions. Arrows 
indicate the positions of the relevant A ze ro points. 



nique could provide much-needed data to constrain the 
BBR shift for optical lattice clocks. 

Measuring light shifts near a A zoro to constrain matrix 
elements has a number of advantages over the approaches 
described earlier, ft is insensitive to absolute calibration 
of the applied optical field, requiring only a stable inten- 
sity. Scalar light shifts are present for all atomic states, 
and A zoro points are found near every atomic transition. 
That proximity increases sensitivity to the matrix ele- 
ment of the nearby transition in addition to decreasing 
sensitivity to uncertainty in hard-to-calculate theoretical 
values. Most importantly, ground (or metastable) state 
light shifts depend only on matrix elements directly to 
higher excited states, not on the coupling of the excited 
states to other states. Light shift measurements are thus 
independent of the branching ratios that limit the life- 
time approach. 

We used this technique to measure the transition ma- 
trix elements between the ground 5si/2 and 6pi/2, &P3/2 
states in 87 Rb by measuring the light shift as a function 
of A near the two A zor o points neighboring the 6p states 
(see the inset in Fig. Q]). We applied a sequence of stand- 
ing wave (optical lattice) pulses to a Bose-Einstein con- 
densate (BEC) containing ~ 3.5 x 10 4 |F = l,mp = —1) 
atoms with no discernible thermal fraction, produced in 
a hybrid optical and magnetic trap similar to [25] . At the 
end of the pulse sequence, we suddenly turn off the trap 
and measure the atom population in each diffracted order 
from an absorption image after 40 ms time of flight. By 
intentionally using a small condensate, we can neglect 
atom-atom interactions, and limit the optical depth so 
that all orders in the image are unsaturated. The mea- 
sured population fractions, normalized to the total num- 
ber in each image, are insensitive to small fluctuations in 
total atom number. 

The lattice light is provided by a frequency-doubled 



diode laser, and is tunable between 419 nm and 424 nm. 
We measure the wavelength to 50 fm (90 MHz) accuracy 
with a wavemeter, calibrated to the known 5s-6p transi- 
tion frequencies [26[ • We form the lattice by reimaging a 
retroreflected beam on the atoms with an incident power 
of 60-120 mW and a waist of ~ 110 //m. The lattice beam 
intensity varies only a few percent across the BEC, whose 
Thomas-Fermi radius is 18 /im transverse to the beam. 
Lattice alignment is performed at 421.700 nm where the 
light shift is reasonably large, ~ 7 Er (where the re- 
coil energy Er = h 2 k 2 /(2M) for wavevector k — 2n/\, 
lattice laser wavelength A, and atomic mass M). The 
input beam is carefully aligned to the BEC (with the 
retroreflected beam blocked) by minimizing the trans- 
verse displacement of the BEC induced by the beam. 
The retroreflected beam is then aligned by maximizing 
diffraction efficiency. Lattice depth measurements at the 
5% uncertainty level are consistent with no drift during 
a data set, and the alignment is typically stable from day 
to day. 

To minimize wavelength-dependent steering of the 
laser, we couple the light through an optical fiber and 
focus the output of the fiber onto the BEC. We mea- 
sured the residual steering after the fiber by monitoring 
the position of the beam 1 m downstream from the fiber 
tip. The observed deflection of the beam was < 30 /xrad, 
which corresponds to at worst a 1% decrease in average 
beam intensity at the BEC and is included in our un- 
certainty analysis. Additionally, we monitor the power 
in each pulse sequence with a photodiode to account for 
different laser powers. 

Lattice polarization also affects the measured light 
shift. Linearly polarized light has a purely scalar light 
shift, but elliptical polarization will have both vector and 
scalar light shifts for atoms with \mp\ > (2?]], shifting 
the position of A zcro - Although we used a Glan-laser po- 
larizer to establish a clean linear polarization, analyzing 
the light after the chamber showed that the vacuum win- 
dows' slight birefringence created ~ 1% cllipticity. To 
account for the possible systematic shift of A zoro , we mea- 
sured the light shift for two orthogonal, linear input po- 
larizations, defined as S and P with respect to the plane 
of the last mirror before the chamber. By symmetry, the 
induced ellipticity from the window is opposite for the 
two input polarizations. We calculated that averaging 
the S and P measurements cancels the contribution of 
the vector light shift to a negligible 10 fm shift in A zcr o- 

Optical lattice depths (light shifts) are commonly mea- 
sured by observing the population oscillations between 
momenta 2lhk, for integer I, as a function of lattice 
pulse duration (outside the Raman-Nath regime). For 
lattice depths below 1 Er, the diffracted populations are 
small, and the oscillation period saturates to h/(AER). 
At 0.05 Er, less than 0.008% of atoms are diffracted for 
a single lattice pulse. In order to increase the signal, 
we employ a multi-pulse diffraction sequence that coher- 
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FIG. 2. A Bloch sphere representation of the state vector's 
evolution during a pulse sequence. The state vector starts in 
\0hk) (north pole). Lattice (free) evolution is shown in blue 
(red). Alternating lattice and free evolutions with duration 
h/(8En), the diffracted population is maximized as described 
in the text. This sequence provides 50% population in \±2hk) 
after four pulses, corresponding to a lattice depth of ~ 1-Er, 
much larger than depths we actually measured. 



ently adds the effect of each pulse, similar to quantum 
resonances in delta-kicked rotors H, 
tices where only orders with 



< 1 



29J. For weak lat- 
are populated, a 

Bloch sphere picture provides intuition [3(3] (Fig. O. For 
initially stationary atoms, the \l\ = 1 orders are equiva- 
lent and can be represented as a single state, which along 
with the I = state are the two poles of the Bloch sphere. 
Free evolution corresponds to precession about the verti- 
cal axis with period h/ (4Pr), called the Talbot time [3ll | . 
During a weak lattice pulse, the precession axis is tilted 
by an angle proportional to the lattice depth with essen- 
tially the same h/{AEn) period. Alternating lattice and 
free evolution with duration h/(8Ef/,) (one half preces- 
sion period, ~ 9.75 /is at 423 nm) efficiently increases 
the diffracted population (see Fig. [2J. 

To accurately determine the lattice depths from the 
measured populations, we numerically diagonalized the 
lattice Hamiltonian in the plane wave basis including or- 
ders up to \l\ =3. From this we calculate the first-order 
diffracted population fraction Pi as a function of pulse 
number n p and lattice depth Vb in units of Er. In the 
weakly diffracting limit (n p Vo <C 4), Pi oc (n p Vo) 2 , show- 
ing that multiple weak pulses are equivalent to one pulse 
n v times stronger. Inverting the exact numeric result for 
Pi, we extract the lattice depth Vo(n p , Pi) from the mea- 
sured population fraction Pi = (N-i + Ni)/N to t, with Ni 
the number of atoms in the Ith. order. Fig. [3] shows Pi 
versus n p at fixed Vq. The agreement with theory is ex- 
cellent for n p up to 18. We attribute deviations at larger 
n p to the reduced overlap of the diffracted components, 
which travel a significant fraction of the size of the BEC, 
during the pulse sequence. 

To determine the polarizability, we extract the raw lat- 
tice depth Vo(n p ,Pi) for each image, which is then di- 
vided by the measured photodiode power for that pulse 
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FIG. 3. Diffracted population at fixed laser power and wave- 
length. As we vary the number of pulses, the diffracted pop- 
ulation (our signal) initially increases as rip. The solid line is 
is the theoretical calculation with only Vb as a fit parameter. 
For accurate light shift measurements, we only used ri v < 15. 



sequence. Fig. 2] shows the extracted polarizability as a 
function of A for four data sets: two orthogonal polar- 
izations, S and P, near each A zor o- The lattice is aligned 
before each data set, and measurements after confirm 
consistent alignment within a data set. The measure- 
ment is repeated up to 20 times at each wavelength [32| , 
and the error bars in Fig. [4] represent purely statistical 
uncertainty of la. For both sets around the 421 nm zero, 
we used 15 pulses, and for both sets around 423 nm we 
used 8 or 11 pulses. 

We fit the polarizability to the expression calculated 
in (H: 
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where the d npj = (npj\\D\\5s) are reduced dipole matrix 
elements and u> is the lattice laser frequency. The matrix 
elements de Pl/2 and de P3/2 are fit parameters, while ds Pl/2 
and d$ P3/2 are fixed by the experimental lifetime mea- 
surement [33|]. The 7p and 8p matrix elements and the 
contribution C C orc, which includes the polarizability of 
the core electrons and their interaction with the valence 
electron, are calculated in [22|. We used an improved 
calculation of Ctaii, the valence contribution for n > 8, 
with reduced uncertainty (see supplemental material). 

Fitting the data to Eq. [JJ around a given A zcro con- 
strains the relationship among the parameters in the 
equation. For example, the position of the A zero near 
421 nm, which is between the two 6p lines, depends 
strongly on the ratio of the 6p matrix elements, Rq p = 
^6p 3/2 /^6pi /2 j but more weakly on their absolute value. 
On the other hand, the position of the A zcro near 423 nm, 
which is red-detuned of both 6p lines, depends more on 
the average value of the 6p elements (in relation to the 
5p elements). Fits around a single A zcro constrain 6p 
elements with respect to the non-6p parameters (d np , 
Ctaii and C corc ), and the extracted values are sensitive to 
their uncertainties. Including data around multiple A zero 
points further constrains the fits. For our case, where the 
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FIG. 4. Measured polarizability (arb. units) versus wave- 
length for two linear, orthogonal input polarizations S (red) 
and P (blue). Each point is an average of up to 20 mea- 
surements, shown with la statistical error bars. The solid 
lines (nearly indistinguishable) are fits to the expected form 
(Eq. [TJ with the 6p matrix elements as free parameters; re- 
duced x 2 f° r the S and P fits is 4.33 and 15.7, respectively. 
Additionally, we allow a separate amplitude about each zero 
to account for different laser power. Fit residuals are shown 
at the top. 

two A zcro 's are close to each other, the effect of the non- 
6p contributions is essentially identical for the two zeros, 
so that simultaneous fits to both zeros accurately deter- 
mines i?6pj independent of the other parameters. With 
the addition of sufficiently accurate theoretical and ex- 
perimental values for the other d np , Ctaii and C coro , we 
can also accurately determine absolute values for de Pl , 2 
and d 6p3/2 . 

We simultaneously fit the 421 nm and 423 nm data 
sets for a given polarization. Due to possible differences 
in alignment and total power between sets, we fit with 
an independent amplitude, A\, around each A zcro giving 
a total of four fit parameters: A421, ^-A23,d6 Pl/21 de P3/2 . 
The fit values are then averaged for the two polariza- 
tions. (The small difference in extracted matrix elements 
for S and P polarization is consistent with the 1% ellip- 
ticity discussed earlier.) Table Q] summarizes the contri- 
butions to our final uncertainty in our measured values 
for i?6 P , de Pl/2 and de P3/2 . The statistical uncertainty 
in Rq p , de Pl/2 and de P3/2 are roughly the same, of order 
0.1%. As expected, the contribution of the uncertainties 
of the other fixed parameters in Eq. Q] to R 6p is negligi- 
ble, providing a good comparison with theory, which can 
predict ratios more accurately than matrix elements (see 
supplemental material). The uncertainties in de Pl/2 and 



TABLE I. Absolute uncertainty contributions (in eao X 1CP 4 ) 
for the 5s-6p matrix elements and their ratio (xlCP 4 ). Note 
the insensitivity of Rg p to uncertainty in the fit parameters. 
Total uncertainty is summed in quadrature. Additionally, our 
5s-6p matrix elements are compared to the theoretical values 
(in ea ). 



Contribution 8de Pl/2 8d$ P3/2 8Rq p 

statistical 1.60 1.71 9.99 

d 5pi/2 0.84 1.35 0.005 

d 7pi/2 0.08 0.13 0.011 

d Spi/2 0.02 0.04 0.002 

npi/ 2 tail 0.56 0.92 0.029 

d 5p3/2 1.77 2.87 0.007 

d 7p3/2 0.22 0.36 0.031 

d Sp3/2 0.06 0.10 0.006 

np 3/2 tail 2.01 3.28 0.104 

core 1.25 2.05 0.065 

alignment drift 8.00 6.31 20.4 

Total 8.74 8.29 22.7 



theoretical value [22] 


0.325(9) 


0.528(13) 


1.624(7) 


our results 


0.3236(9) 


0.5230(8) 


1.616(2) 



de P3/2 do depend on the other parameters, and the two 
largest contributions are from the theoretical uncertainty 
in the np 3 / 2 component of Ctaii and the experimental un- 
certainty in d§ P3/2 ; these contribute at the 0.06% level. 
To account for potential lattice alignment drift, we sim- 
ulated the effect of a 5% linear drift in the extracted 
polarizability across a data set, which resulted in an ad- 
ditional 0.2% uncertainty in d 6pi/2 , and 0.1% in R 6p and 
d 6 p 3/2 \M- 

Our light shift cancellation measurement technique ac- 
curately determines the ratio of the 6p matrix elements, 
with Rq p = 1.616(2). The matrix elements are de Pl/2 — 
0.3236(9) and de P3/2 = 0.5230(8), in excellent agreement 
with the theoretical values (Table |T| . From these we de- 
termine values for the A zero points: 421.075(3) nm and 
423.018(8) nm. 

We have presented the first experimental measurement 
of the 5s-6p dipole transition matrix elements in Rb. Our 
technique of measuring the relative light shift near A zcro 
is applicable to many important atoms. At a minimum, 
such measurements constrain matrix elements and pro- 
vide benchmarks needed for atomic theory. Absolute ma- 
trix element values can be bootstrapped from a single, 
well-known matrix element. Additionally, applying our 
measurement technique to metastable states, such as the 
3 Po clock state of Yb or Sr, could provide information 
about excited-excited state matrix elements. 

We acknowledge useful discussions with G.K. Camp- 
bell, L.A. Orozco, and W.D. Phillips. This work was 
partially supported by the nsf-phy 1104472 and the ARO 
with funding from DARPA's OLE program and ARO 
Grant W911NF0910216 (DURIP). The work of MSS was 
performed under the sponsorship of the US Department 



5 



of Commerce, NIST. 



* cherold@umd.edu 
K. Beloy, J. A. Sherman, N. D. Lemke, N. Hinkley, 
C. W. Oates, and A. D. Ludlow, |arXiv e-prints (2012)| 
arXiv:1208.0552vl [physics. atom-ph] 
J. A. Sherman, N. D. Lemke, N. Hinkley, M. Pizzo- 
caro, R. W. Fox, A. D. Ludlow, and C. W. Oates, 
|Phys. Rev. Lett. 108, 153002 (2012)] 
S. G. Porsev, A. D. Ludlow, M. M. Boyd, and J. Ye, 
Phys. Rev. A 78, 032508 (2008). 

S. G. Porsev, K. Beloy, and A. Derevianko, 

|Phys. Rev. Lett. 102, 181601 (2009)] 

A. A. Vasilyev, I. M. Savukov, M. S. Safronova, and 

H. G. Berry, |Phys. Rev. A 66, 020101 (2002) | 

M. S. Safronova, U. I. Safronova, and C. W. Clark, arXiv 

e-prints (2012), [arXiv:1206.7115 [physics.atom-ph]| 

A. V. Gorshkov, A. M. Rey, A. J. Daley, M. M. Boyd, 

J. Ye, P. Zoller, and M. D. Lukin, Phys. Rev. Lett. 102, 

110503 (2009). 

M. Saffman and T. G. Walker, 



Phys. Rev. A 72, 02 2347 (2005)] 

H. P. Summers, N. R. Badnell, M. G. O'Mullane, A. D. 
Whiteford, R. Bingham, B. J. Kellett, J. Lang, K. H. 
Behringer, U. Fantz, K.-D. Zastrow, S. D. Loch, M. S. 
Pindzola, D. C. Griffin, and C. P. Ballance, Plasma Phys. 
Control. Fusion 44, B323 (2002). 

D. A. Verner, E. M. Verner, and G. J. Ferland, 
|At. Data Nucl. Data Tabl. 64, 1 (1996) 



J. C. Pickering, R. Blackwell- Whitehead, A. P. 
Thorne, M. Ruffoni, and C. E. Holmes, 

|Can. J. Phys. 89, 387 (2011)| 

M. S. Safronova, M. G. Kozlov, and C. W. Clark, 
|Phys. Rev. Lett. 107, 143006 (2011)] 
S. G. Porsev and~ A. 
|Phys Rev. A 74 020502 (2006)| 
K. Aikawa, AT Frisch, M. Mark, 
A. Rietzler, R. Grimm, and F 



Phys. Rev. Lett. 108, 210401 (2012) 
M. Lu, N. Q. Burdick, a m i L> . 
|Phys. Rev. Lett. 108, 215301 (2012) 



M. 



Saffn 



and 



K. 



Derevianko, 

S. Baier, 
Ferlaino, 

Lev, 

M0lmer, 

O. Dulieu, 



|Phys. Rev. A 78, 12336 (2008) 
N. Bouloufa, AT Crubellier, and 
|Physica Scripta T134, 014014 (2009)] 
E. Gomez, S. Aubin, L. A. Orozco, and G. D. Sprouse, 
|J. Opt. Soc. Am. B 21, 2058 (2004)| 
J. A. Sherman, T. W. Koerber, A. Markho- 
tok, W. Nagourney, and E. N. Fortson, 

|Phys. Rev. Lett. 94, 243001 (2005)] 
J. A. Sherman, A. Andalkar, W. Nagourney, and E. N. 
Fortson, |Phys. Rev. A 78, 052514 (2008) | 
B. K. Sahoo, L. W. Wansbeek, K. Jungmann, and 
R. G. E. Timmermans, |Phys. Rev. A 79, 052512 (2009)| 
B. Arora, M. S. Safronova, and C. W. Clark, 
|Phys. Rev. A 84, 043401 (2011)| 

These are the "tune-out" wavelengths proposed for 
species-selective traps by L. J. LeBlanc and J. H. Thy- 
wissen, [Phys. Rev. A 75, 053612 (2007)| 
J. Ye, H. J. Kimble, and H. Katori, 



|Science 320, 1734 (2008) | 
[25] Y.-J. Lin, A. R. Perry, R. L. Compton, I. B. Spi elman, 

and J. V. Porto, |Phys. Rev. A 79, 063631 (2009) | 
[26] A. Shiner, A. Madej, P. Dube, and J. Bernard, 
|App. Phys. B 89, 595 (2007)] 

I. H. Deutsch and P. S. Jessen, Opt. Commun. 283, 681 
(2010). 

F. M. Izrailev and D. L. Shepelyanskii, 
|Theor. Math. Phys. (Engl Trans.) 43, 553 (1980) | 

Robinson, C. F. 
and M. G. Raizen, 



[27] 
[28] 
[29] 

[30] 
[31] 



F. L. Moore, J. 
Bharucha, B. Sundaram, 
Phys. Rev. Lett. 75 4598 (1995)] 

S. Wu, Y.-J. Wang, Q. Diot, and M. Prentiss, Phys. Rev. 
A 71, 043602 (2005). 

L. Deng, E. W. Hagley, J. Denschlag, J. E. Simsarian, 
M. Edwards, C. W. Clark, K. Helmerson, S. L. Rolsto n, 
and W. D. Phillips, [Phys. Rev. Lett 83, 5407 (1999)] 
[32] Some points exhibited decreased diffraction, correlated 
with the condensate position in each image, due to resid- 
ual condensate motion. Thus we excluded the 22% of im- 
ages where the condensate was more than 15 /im from the 
average position for each data set. The extracted matrix 
element is insensitive to the exclusion value. 

[33] XL Volz and H. Schmoranzer, 

|Physica Scripta T65, 48 (1996) | 
[34] That reduced \ z f° r the fits in Fig. [3] is larger than 1 is 
accounted for by the additional uncertainty in alignment 
drift; rescaling the statistical uncertainty to give reduced 
\ 2 of 1 produces similar error bars. 
[35] M. S. Safronova, W. R. Johnson, and A. Derevianko, 

Phys. Rev. A 60, 4476 (1999). 
[36] Y. Ralchenko, F. C. Jou, D. E. . Kelleher, . E. Kramida, 
A. Musgrove, J. Reader, W. L. Wiese, and . Olsen, "Nist 
atomic spectra database," (2005), (version 3.1.2). [On- 
line]. Available: http://physics.nist.gov/asd3 [2007, Au- 
gust 29]. National Institute of Standards and Technology, 
Gaithersburg, MD. 
[37] J. Sansonetti, W. Martin, and S. Young, "Handbook of 
basic atomic spectroscopic data," (2005), (version 1.1.2). 
[Online] Available: http://physics.nist.gov/Handbook 
[2007, August 29]. National Institute of Standards and 
Technology, Gaithersburg, MD. 



6 



Supplemental Material 

Calculation of the magic-zero (also called "tune-out") wavelengths where the frequency-dependent polarizability 
vanishes was discussed in detail in Ref. 22]. Briefly, the frequency-dependent scalar polarizability, ao(u>), of an 
alkali- metal atom in its ground state v may be separated into a contribution from the core electrons, a C orc, a core 
modification due to the valence electron, a vc , and a contribution from the valence electron, q;q(w). The valence 
contribution to the static ac polarizability is calculated using the sum-over-states approach: 

„, x 2 ^ (k\\D\\v) 2 {E k ~E v ) 

a ° M = ^TT)V (E k - Ev r-^ ' (2) 

where (k \\D\\ v) is the reduced electric-dipole (El) matrix element. Because of the rapid convergence of the sum over 
intermediate k states in Eq. @, we separate the valence state polarizability into two parts, a ma i n , containing the 
contributions from the few lowest np states, and the remainder, ataii- The experimental data [33j are used for 5s — bp 
matrix elements. The 5s — 6p matrix elements are determined in the present work. For all other main terms, we use 



the relativistic all-order values 35] of the matrix elements and the experimental values of the energies [36. 37]. In 
Ref. [22| , the contributions with n = 5 — 8 were calculated accurately and tail was calculated in the Dirac-Hartree-Fock 
approximations. To decrease uncertainty in the tail, we calculate n — 9 — 12 contributions using all-order matrix 
elements and experimental energies. The contributions to the frequency-dependent polarizability of the ground state 
of Rb at A zcro = 421.07476 nm and A zero = 423.01787 nm magic-zero wavelengths are given in Table [TT1 We list a corc 
and a vc contributions together in the row "Core" . The uncertainties in the values of the magic-zero wavelengths are 
taken to be the maximum difference between the central value and the crossing of the ao ± Sao with zero, where 
Sao is the uncertainty in the ground state polarizability value at that wavelength. The resulting values of A zor o 
421.075(3) nm and 423.018(8) nm are in excellent agreement with 421.08(3) nm and 423.05(8) nm predictions of 
Ref. [22] ]. respectively. 

The lowest order and all-order values of the ratio R of the 5s — 67)3/2 and 5s — 6pi/2 matrix elements are listed in 
Table [TTT1 In the relativistic all-order method, all single-double (SD) or single-double and partial valence triple (SDpT) 



excitations of the Dirac-Fock (DF) wave function are included to all orders of perturbation theory [35j . Additional 



semi-empirical scaling of the all-order values was carried out for both SD and SDpT values. The scaled SD values 



TABLE II. Contributions to the frequency-dependent polarizability of the ground state of Rb at A ze ro = 421.07476 nm and 
A Z oro = 423.01787 nm. Absolute values of electric dipole matrix elements are expressed in a.u. (eao), and the corresponding 
energy differences are expressed in conventional wavenumber units (cm -1 ). 



Contribution 


\{nps\\D\\5a)\ 




Q?421. 07476 


C*423. 01787 


5f>l/2 


4.231(3) 


12578.95 


-40.599(60) 


-41.122(60) 


6P1/2 


0.3235(11) 


23715.08 


-113.608(771) 


50.840(346) 


7P1/2 


0.115(3) 


27835.02 


0.128(6) 


0.125(6) 


8P1/2 


0.060(2) 


29835.00 


0.024(2) 


0.023(2) 


9f>l/2 


0.037(3) 


30958.91 


0.0080(11) 


0.0079(11) 


IOP1/2 


0.026(2) 


31653.85 


0.0036(5) 


0.0036(5) 


HPl/2 


0.020(1) 


32113.55 


0.0020(3) 


0.0019(3) 


12pi/ 2 


0.016(1) 


32433.50 


0.0012(2) 


0.0012(2) 


(n > 12)pi/ 2 tail 






0.039(39) 


0.038(38) 


Total (n > 8)pi/ 2 tail 






0.054(41) 


0.053(41) 


5P3/2 


5.978(5) 


12816.55 


-83.824(126) 


-84.937(129) 


6P3/2 


0.5230(11) 


23792.59 


228.446(962) 


65.637(276) 


7p 3 /2 


0.202(4) 


27870.11 


0.392(17) 


0.383(16) 


8P3/2 


0.111(3) 


29853.80 


0.082(5) 


0.081(4) 


9P3/2 


0.073(5) 


30970.19 


0.0303(44) 


0.0299(43) 


10p 3 /2 


0.053(4) 


31661.16 


0.0146(21) 


0.0144(21) 


11P3/2 


0.040(3) 


32118.52 


0.0082(12) 


0.0081(12) 


12P3/2 


0.033(2) 


32437.04 


0.0053(8) 


0.0052(8) 


(n > 12)p 3 /2 






0.138(138) 


0.136(136) 


Total (n > 8)p 3 / 2 tail 






0.196(146) 


0.194(144) 


Core 






8.709(93) 


8.709(93) 


Total a 






0.000 


0.000 
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TABLE III. The values of the 5s — 6pj reduced electric-dipole matrix elements in atomic units and their ratio 7?. The first-order 
(DHF) and all-order SD and SDpT values are listed; the label "sc" indicates the scaled values. 







DHF 


SD 


SDpT 


SD SC 


SDpT sc 


Final 


5s 


- G'Pl/2 


0.3825 


0.3335 


0.3300 


0.3248 


0.3338 


0.325(9) 


5s 


~ 6p 3/2 


0.6055 


0.5409 


0.5354 


0.5276 


0.5401 


0.528(13) 


R 


1.5831 


1.6217 


1.6223 


1.6245 


1.6179 


1.624(7) 



(SD SC ) is taken as final. The uncertainty in the values of the ratio is estimated as the the maximum difference of 
the SD, SDpT, and SDpT sc values from the final result. The same procedure was used to determine theoretical 
uncertainty in the values of both matrix elements. 



